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A model of pulsed inductive plasma thrusters consisting of a set of coupled circuit equa- 
tions and a one-dimensional momentum equation has been used to study the effects of 
adding a second, parallel capacitor into the system. The equations were nondimensional- 
ized, permitting the recovery of several already-known scaling parameters and leading to 
the identification of a parameter that is unique to the particular topology studied. The 
current rise rate through the inductive acceleration coil was used as a proxy measure- 
ment of the effectiveness of inductive propellant ionization since higher rise rates produce 
stronger, potentially better ionizing electric fields at the coil face. Contour plots repre- 
senting thruster performance (exhaust velocity and efficiency) and current rise rate in the 
coil were generated numerically as a function of the scaling parameters. The analysis re- 
veals that when the value of the second capcitor is much less than the first capacitor, the 
performance of the two-capacitor system approaches that of the single-capacitor system. 

In addition, as the second capacitor is decreased in value the current rise rate can grow to 
be twice as great as the rise rate attained in the single capacitor case. 

I. Introduction 

P ulsed inductive plasma accelerators are electrodeless space propulsion devices where a capacitor is 
charged to an initial voltage and then discharged as a high-current pulse through a coil, inductively 
coupling energy into the propellant. The field produced by this pulse ionizes the propellant, producing a 
plasma near the face of the coil. Once a plasma is formed it can be accelerated and expelled at a high exhaust 
velocity by the Lorentz force arising from the interaction of an induced plasma current and the magnetic 
field. While there are many coil geometries that can be employed to inductively accelerate a plasma, in 
this paper the discussion is limited to planar geometries where the coil takes the shape of a flat spiral. A 
recent review of the developmental history of planar-geometry pulsed inductive thrusters can be found in 
Ref. [1]. Two concepts that have employed this geometry are the Pulsed Inductive Thruster (PIT) 2,3 and 
the Faraday Accelerator with Radio- frequency Assisted Discharge (FARAD). 4 

Pulsed inductive plasma accelerators possess many demonstrated and potential benefits, 1,3 providing 
motivation for continued investigation. The electrodeless nature of these thrusters eliminates electrode 
erosion, mitigating the lifetime and contamination issues associated with conventional electric thrusters. 
Also, a wider variety of propellants is available for use when compatibility with metallic electrodes is no 
longer an issue. Pulsed inductive accelerators have demonstrated operation on propellants like ammonia, 
hydrazine, and CO 2 , and there is no fundamental reason why they would not operate on other propellants 
like H 2 O. It is well known that for a given propellant pulsed accelerators can maintain constant specific 
impulse J sp and thrust efficiency rjt over a wide range of input power levels by adjusting the pulse rate to 
maintain a constant discharge energy per pulse. In addition, these thrusters have demonstrated operation 
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Figure 1. A) General lumped-element circuit model and B) equivalent circuit model of a pulsed inductive 
accelerator (after Ref. [2]). C) Equivalent electrical circuit model of an accelerator with a second, parallel 
capacitor. 


in a regime where rjt is relatively constant over a wide range of J sp . Finally, thrusters in this class have 
operated at high energy per pulse, and by increasing the pulse rate they offer the potential to process very 
high levels of power to provide relatively high thrust using a single thruster. 

Pulse circuits for inductive thrusters have in the past typically been limited to a simple, ringing resistive- 
inductive-capacitive (RLC) configuration like that shown in Fig. 1 A and redrawn in Fig. IB as an equivalent 
circuit that is easier to model. However, as the field has developed proposed circuit topologies are becoming 
more complex . 5 In this paper, we proceed with an investigation of the circuit shown in Fig. 1 C where a 
second capacitor with value less than or equal to C\ is inserted downstream of the switch. This study was 
motivated by a data set where the efficiency of a pulsed inductive thruster increased when the capacitor C2 
was added . 6 Unfortunately, the value of C\ was also increased when C2 was added and previous work has 
shown that increasing C\ can also increase the efficiency . 7 The authors of Ref. [ 6 ] additionally noted that the 
voltage across C2 could be approximately double the charge voltage across C\ when C2 <C C\. This result 
was interesting because it implied that the voltage and commensurate current rise rate in the coil could be 
increased by adding C2. A higher current rise rate can, in turn, produce stronger electric fields at the coil 
face and potentially lead to better inductive ionization of the propellant . 5 

There exists a 1 -D pulsed inductive acceleration model that employs a set of circuit equations coupled 
to a one-dimensional momentum equation . 2,3 The model has since been nondimensionalized and used by 
Polzin et al 7,8 to define a set of scaling parameters and gain general insight into their effects on thruster 
performance. In this paper we modify the acceleration model to account for the presence of C2 in the system 
and then nondimensionalize the equation set to identify any new nondimensional scaling parameters that 
arise for the new circuit topology. The current rise rate through the coil is computed for various cases, and it 
is used as a proxy measurement to gauge the ability of the coil to inductively ionize the propellant. Finally, 
we examine whether the addition of C2 imposes potential benefits or detriments on thruster efficiency and 
Isp • 


II. Governing Equations 


A. Physical Model 

A circuit-based model of a pulsed inductive accelerator previously developed by Lovberg and Dailey 2 is 
modified to account for the addition of the second capacitor as shown in Fig. 1 C. In the lumped-element 
circuit model shown in Fig. 1 A,B, the external circuit (left side of the figure) possesses a capacitor bank with 
capacitance Ci, external inductance Lo, resistance R e , and acceleration coil inductance Lc . The plasma 
also has an inductance equal to Lc and a resistance R p . The two circuits are inductively coupled through 
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the coil, which acts as a transformer with mutual inductance M . The value of M is a function of the current 
sheet position z. In Fig. 1C, a second capacitor C 2 is added between the primary capacitor bank and the 
inductive coil. This has the effect of splitting the Lq and R e , with the value on the left and right sides of C 2 
linearly scaled by the factor n < 1 and (1 — n), respectively. Physically, we are assuming that the external 
inductance and resistance in each leg is scaled by the position of C 2 relative to C\ and the coil, with n 
approaching unity as C 2 is moved nearer to the coil. 

A set of circuit equations for Fig. 1C are written through the application of Kirchoff’s law to each loop. 
Rearranging the equations and adding a statement for the time-rate of change of charge on each capacitor 
yields the following coupled set of first-order ordinary differential equations: 


dh 

dt 
<U 2 
dt 

dh 

dt 
dV 1 
dt 
dV 2 
dt 


V 1 -V 2 - nR e h 
nL 0 

{ L c V 2 + ( Lch + Mh) ^ - (1 - n)R e Lch - R P MI 3 j / {L c [(1 - n) L 0 + L c ] - M 2 } , 


T dM dl 2 T 
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_h_ 

Ci 

h-h 

c 2 


where V\ and V 2 are the voltages on capacitors C\ and C 2 , respectively. 

The inductance of a planar, spiral inductive coil coupled to a plasma current sheet is difficult to compute 
analytically. However, based on experimental measurements it has been found 2 that the mutual inductance 
can be modeled using the simple function 


M = Lcex p(~), (2) 

which, when differentiated results in the following equation that governs the changing mutual inductance of 
the circuit: 

dM Lc ( z \ dz 

dt 2zo ^ \ 2zo J dt 

As the current sheet moves forward, it entrains and accelerates any enocuntered gas. The propellant 
mass in the current sheet as a function of time can be written as 

m(t) = mo + / pav z dt , (4) 

Jt = 0 

where pA = pA{z(t)) is the linear mass density distribution and v z is the sheet velocity. The term mo 
represents the initial mass of propellant in the sheet while the integral term represents the mass accumulated 
by the sheet as it moves away from the acceleration coil. 

The momentum equation for this system can be written as 

tMl. exp {-z/z 0 ) = p A v 2 z +m(t ) ^. (5) 

The left hand side represents the self-field electromagnetic force generated through the interaction of the 
current and the magnetic field. The first term on the right hand side in Eq. (5) is the momentum investment 
associated with ‘snowplowing’ propellant (i.e. accelerating the newly entrained propellant encountered by 
the sheet from rest to the sheet speed) while the second term involves further acceleration of the already 
entrained propellant. 

The model lacks a plasma model to calculate temperature and resistivity and it does not contain any 
plasma ionization, or ‘breakdown’, model, instead simply assuming that there is plasma present at t = 0. 
These shortcomings were previously discussed in more detail in Ref. [7]. As in that previous work, the purpose 
here is to use the model to search for relevant nondimensional scaling parameters and gauge their general 
effect on performance, instead of looking for an exact matching of experimental data to the performance 
model. However, we do note that even with the shortcomings, the 1-D, circuit-based acceleration model has 
shown good qualitative and quantitative agreement with experimental data. 1,3 
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B. Nondimensionalization of the Model 

The nondimensionalization strategy previously employed with success in Ref. [7] can also be used in this 
case. The following dimensionless terms are selected to aid in the conversion of the dimensional equations 
to a nondimensionalized set: 


'i = vM h 


1 Lo 
V 0 \ C 

t 


n = -\i^h 

t* = 


(6) 


V LoC 

Other dimensionless terms that arise in the course of nondimensionalizing the equation set are: 


2 * y* _ Yl 

Z ~ > ’'I T/ > 

Zo yQ 

M * VL^C 

M* = —, v* = v z . 

L c z o 


T „ v 2 
^ = vv 


(7) 


The current sheet mass can be nondimensionalized as 

r 1 

/o 


m* =ml + [ p*f(z *) v* dt* 

Jo 


(8) 


where m g = mo /mbit, p* = po^o/^bit? and mbit is the total propellant mass per pulse. 

Using these nondimensional terms, equations (1), (3), (4), and (5) are written in dimensionless form as, 
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conditions for solving the nondimensional set of equations are as follows. 


(9a) 

/ \(l-n)L* + 1) - (M*) 2 ](9b) 

(9c) 

(9d) 

(9e) 

(9f) 

(9g) 

(9h) 
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1*1 (0) = 0, 

M* (0) = 1, 

I 2 (0) = 0 , 

(0) = 0 , 
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C. Scaling Parameters 

During the course of nondimensionalization, several new terms appear in Eqs. (9). These new terms are the 
relevant scaling parameters of the system of equations and are defined as 



Vh = Re 




cl yjLc 

2 mbit zl ’ 


( 11 ) 


which are similar to the parameters previously found for the single capacitor case. 7 We proceed with a brief 
discussion of the meanings of the scaling parameters. 


1. Inductance Ratio L* 

A pulsed inductive accelerator circuit possesses an external inductance Lq and an acceleration coil inductance 
Lc • During a current pulse the moving plasma increases the circuit’s inductance from Lq to LqPLq (he., 
Lc'=AL). The fractional change of inductance, (L*)~ 1 =AL/ Lq, h 1 a pulsed electromagnetic accelerator 
provides a measure of efficiency, 9 as this ratio is indicative of the fraction of energy that can be deposited 
into electromagnetic acceleration of the gas. In an efficient pulsed inductive accelerator the value of L* must 
be much less than unity or Lc Lq. 


2. Dynamic Impedance Parameter a 

The dynamic impedance parameter a can be recast as the product of several important ratios: 

ClVgLc 1 CW 0 2 /2 / 2WZ^ \ 2 

2mbit^o 87r 2 m b it^/2 \ L 0 /L ) 

where L is the dynamic impedance, defined as v z L', and L' is the effective inductance per unit length, 
equal to Lc/zq. The ratio of the initial stored energy (C\Vq / 2) to the plasma kinetic energy (rabit^/2) 
is recognized as the inverse of thrust efficiency ry x and will always be greater than one. The final term is 
the ratio of the natural period of the driving external circuit, 2tt^/LqCi, to the time interval, Lq/L, over 
which the plasma’s motion increases the inductance of the circuit by one unit of Lq. The former is the 
timescale on which the external circuit naturally operates, while the latter is the timescale on which the 
current sheet remains electromagnetically coupled to the acceleration coil. For a given configuration, there 
exists an optimum value of a where the electromagnetic coupling timescale is matched to the frequency of 
the external circuit, allowing for optimum transfer of stored electrical energy into directed kinetic energy. 7 


3. Circuit Parameters t/q , ^2 and C 


The critical resistance ratios ipi and ^2 were shown in previous work to control the nature of the current 
waveforms. 7 The new parameter C also has an effect on the nature of the current waveforms. Proceeding 
with an analysis similar to that performed in Ref. [7], we search for limiting-case solutions to Eqns. (9a)- (9e). 
Decoupling the current sheet dynamics (i.e., the acceleration and sheet motion) from the problem allows us 
to apply the condition 


M* = 1, 


dM* 

dt* 


= 0, 


which simplify the circuit equations. Under these assumptions, the circuit equations can be combined and 
rewritten as 


d 2 It 

dt* 2 
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If the right side of the third equation in (13) is small, then the induced current in the plasma mirrors the 
current in the coil 


This is the case where the coil and plasma currents are well-coupled, which is expected in efficient accelerators. 
The third equation of (13) can be substituted into the second to yield 


a 


d?i* 2 

n) 


+ [(i 


cfT* 


+ Cl 2 * = CI{ 


■ 4> 2 2 L*I* 


(14) 


This second-order non-linear ODE is similar to the one found in Ref. [7]. If the last term on the right-hand 
side is comparitively small, then the forcing term contains only /*, which will never be small and cannot be 
neglected. We can gain insight by writing the homogeneous solution 


r* 

J 2 


Aq exp (— Tt*) sin 



(15) 


where 

*= (^ i + (i^)) /2 - 

In this case, the damping is a function of the critical resistance ratios while the natural frequency is addi- 
tionally a function of the parameter C. 

The term on the right hand side of the first equation in (13) will also never be small. However, we can 
write the homogeneous solution for I* as 


I* = Bq exp 



sin 


'(C + 1) 



(16) 


Though these are only the homogeneous solutions to the simplified electrical response equations, there 
are several observations that can be made as a result of this exercise. The parameters i/j i, ^ 2 , and C control 
the electrical responses in the circuit. Increasing the parameter C increases the ringing frequency of both 
I* and /|, with the frequencies approaching each other as C becomes very large. To neglect the non-linear 
terms presented by the current J|, the values of L* and ^2 must be such that 

^\L* « 1 , 
ip 2 L* <S 1. 

Finally, Eq. (9e) can be used to show that as C approaches infinity, the derivative dV£ / dt * will only remain 
bounded if I* = . Physically, this is the case where C 2 is removed from the system, permitting the recovery 

of the response for the circuit shown in Fig. IB. 


III. Nondimensional Solutions 


A. Solution Strategy 

The set of coupled first-order ODEs given in Eqs. (9) can be solved numerically once the mass distribution 
and the set of scaling parameters given in Eqs. (11) are specified. The performance metrics chosen for this 
study are the exhaust velocity, u*, and the thrust efficiency, which is written in terms of nondimensional 
parameters as 

to* ( v*. f 

1h = -2L^r- 

Note that this corrects the equation given for efficiency in Ref. [7] where the square term was left off u* 
(though it was included in the analysis contained in the reference). 

In solving any set of first-order (in time) differential equations, it is important to know when the time 
histories of the computed variables (specifically u* in our case) should be queried to calculate performance. 
This question is, in fact, critical to the evaluation of these accelerators. For our non-dimensional model, the 
integration period will end when one of the following two conditions is reached: 
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1. The end of the first half-cycle of the accelerator coil discharge is reached and the current I* reverses 
in sign. 

2. The sheet travels three characteristic lengths, z* = 3. 


The first condition is based on the fact that when the accelerator current goes through zero, it is going 
though a point of high dl/dt. While the acceleration model doesn’t incorporate any ionization physics, it is 
well known that a new current sheet can form at the face of the coil during current reversals, causing what 
is known as a “crowbar discharge” . If this occurs, the initial current sheet ceases to undergo acceleration. If 
at the time of current reversal the final current sheet mass m J is less than the injected propellant mass bit 
ra£ it , conservation of momentum can be used to correct the end-of-calculation final value of velocity (v*)f 
to give the velocity as 


v 


* 


mn< 


mn 


«)/• 


bit 


(18) 


The second condition stems from the existence of a finite axial distance between the current sheet and the 
coil, beyond which the two are essentially decoupled electromagnetically. The separation distance z*= 3 is 
chosen as our cutoff for electromagnetic coupling as it represents an inductance change in the circuit of 95% 
of the coil inductance. Above this cutoff value, the integration yields no significant change to the calculated 
performance. 

For the present study, a triangular (linearly decreasing) mass distribution is used. This is given by 


Pa = 


Po (1 - z*/5* m ) 

0 


z* <5 
z* >5 


* 

m 

* 

m') 


where = 2/^m, = 3m/z o, and S m is the (dimensional) depth of the propellant layer. 


B. Solutions 

1. Electrical Response as a Function of C and n 

We first present representative waveforms showing the electrical behavior of the accelerator as a function of 
C in Fig. 2 and n in Fig. 3. The former gauges the effect of the capacitance ratio C while the latter gives 
insight into the effect of the position of capacitor C 2 in the line relative to C\ and the coil. Waveforms for 
and dPf/dt* are presented, where the latter provides a measure of the inductive voltage across the 
coil and are being used in this study as proxies for the ability of the coil to inductively ionize propellant. 

In Fig. 2 we observe that currents I* and have the standard damped-sinusoidal response of an LRC 
circuit with a higher-frequency waveform superimposed on the response. The modulation in all the waveforms 
increases in frequency and decreases in amplitude as C is increased. For C = 1 we observe a nearly resonant 
energy transfer between C\ and C 2 with the rise rate of the current in the coil relatively low in comparison 
with the other cases. As C is increased to 100, the maximum value of dl^/dt* approaches 2. 

The changing value of n appears to most affect the response of the current I* in Fig. 3, with the amplitude 
of the oscillation superimposed on the wave decreasing with increasing n. We observe that for smaller n 
values dl^/dt* oscillates much longer and, like I{ , with greater amplitude, having the greatest peak value 
for the smallest n displayed. It is easiest to observe in Fig. 3G,H that the overall average behavior of I* and 
I 2 are roughly equivalent ( /* ~ J|) for the cases where C equals 100, which was predicted in the limiting 
case solutions for C sufficiently large. The form and magnitude of the waveform for 1 1 does not appear to 
change much with n, implying that the acceleration process is generally unaffected by the position of C 2 in 
the line. For the remainder of this study we use n = 0.2 to obtain greater current rise rates. 


2. N ondimensional Performance 

Contours of constant computed performance ( 7 ft , v*) and dPf/dt* found through solving the governing 
equations for a range of C and a values are presented in Fig. 4A-C. For a fixed C the efficiency possesses a 
local maximum with respect to <a, which is consistent with the interpretation of a as a dynamic impedance 
matching parameter. The efficiency at a fixed a initially increases with increasing C, eventually asymptoting 
at C between 10 and 100. The velocity increases with a and like the efficiency it also asymptotes as C is 
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f f f 

Figure 2. Waveforms for 7j% 1%, and dPf/dt* showing the variation in each for A)-C) n = 0.2, C = 1, D)-F) 
n = 0.2, C = 10, and G)-I) n = 0.2, C = 100. 









f 


f 


Figure 3. Waveforms for 7J, 7^, and dl^/dt* showing the variation in each for A)-C) n = 0.2, C = 100, D)-F) 
n = 0.5, C = 100, and G)-I) n = 0.8, C = 100. 
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increased. The current rise rate is a weak function of cq but increases very quickly as C is increased until 
asymptoting just below a value of dl | /dt* equal to 2. 

It is worth comparing solutions to the two-capacitor case (Fig. 1C) with solutions to the single-capacitor 
case (Fig. 1A,B) to gauge the effect of the added capacitor C 2 on accelerator performance. Curves of 
computed performance ( 7 ^, v*) and the maximum value of dl^/dt* for C equal to 1.5, 5, 10, and 100 are 
presented in Fig. 4D-F. An additional curve on each graph appears for the single-capacitor case (with the 
maximum value of dl * /dt* given in Fig. 4F). We observe that there is not a large difference in the terminal 
exhaust velocity v* of the single and two-capacitor cases for any value of C. The efficiency graph shows 
that the two-capacitor case asymptotes to the single-capacitor solution as C approaches 100, implying that 
a capacitor C 2 can be added without detrimentally affected thruster performance. Additionally, as C is 
increased the maximum current rise rate through the coil in the two-capacitor case asymptotes to a value 
that is twice the single-capacitor case. 




Figure 4. Contour plots of A) efficiency, B) nondimensional velocity, and C) maximum value of dl^/dt* as a 
function of varying C and a. Values of ^ 1 , ^ 2 ? and L* are 0.05, 0.13, and 0.121, respectively. Comparison of D) 
efficiency, E) nondimensional velocity, and F) maximum dl* /dt* through the inductive coil as a function of a 
for various values of C and for the single-capacitor case. The values of maximum dPf/dt* are given in F) except 
for the single capacitor case, where the maximum dl^/dt* is shown. 
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IV. Discussion 


The data presented in the previous section helped gauge the effect on performance that the addition of 
capacitor C2 would have. Comparing the efficiency for the single and two-capacitor cases demonstrates that 
the performance above a certain value of C is not adversely affected by the the addition of the capacitor 
C2. This result is somewhat expected because the overall behavior of current (seen in Fig. 3 ) is relatively 
well-behaved, with only a small modulation superimposed on top of the main waveform, as opposed to the 
much larger excursions observed in Fig. 2 for smaller values of C. In addition, at C of 100 the waveforms for 
I* and 1 2 are approximately equivalent, implying that the presence of the capacitor C2 is having little effect 
on the overall current through the inductive coil that is inducing fields in the plasma and acting on the gas 
to accelerate it. 

As C is increased, the value of dl^/dt* asymptotes to twice the value of dl* /dt* through the coil in the 
single capacitor case. The reason for this can best be understood by focusing on the first current loop (ii) in 
Fig. 1C. If there is no resistance in the loop, when C\—C2 (C = 1 ), the charge from C\ transfers completely 
to C2, with the current I* peaking when the voltage on each capacitor is equal and going to zero as the 
voltage on C2 equals the original voltage on C \ . Consequently, the voltage across the coil is still roughly the 
same value as it would have been if C2 was not present. This neglects resistive losses and doesn’t account 
for the fact that much of the charge, instead of going through the coil from C2, instead transfers back to C\ 
robbing potential current and dl^/dt* through the coil (a situation illustrated in Figs. 2A-C.). 

We can best understand the case where C\ C2 by further simplifying current loop I * so it consists 
only of Ci, C2, and uLq. (Neglecting the current loops and 1 3 and their associated lumped electrical 
elements does not detract from the reasoning that follows and serves to simplify the explanation.) Initially, 
the voltage on Ci changes very little as the discharge progresses. The maximum current 1 ^, corresponding 
to dl*/dt * = 0, occurs when the voltage on C2 is approximately equal to the initial charge voltage on Ci. 
As time progresses further the voltage on C2 continues to rise, eventually reaching a theoretical maximum 
value that is approximately twice the initial voltage on C\ when the current I* = 0 and the value of dl*/dt* 
is a local minimum. Translating this to the full circuit, we observe that this would imply a voltage across the 
coil that is approximately twice the voltage possible in the single capacitor case, yielding a commensurate 
doubling of dl* /dt* through the coil. 


V. Conclusions 

We have used a nondimensional pulsed inductive acceleration model to study the effects of adding a 
second, parallel capacitor into the circuit topology. Inserting the second capacitor into the system was 
motivated by limited data showing that the voltage appearing across the accelerator coil could be increased 
significantly over the initial charge voltage. Increasing the voltage across the coil has been shown to be 
important in the process of inductively ionizing the propellant. The nondimensionalization of the modified 
model led to identification of new scaling parameters that are unique to the particular topology. The physical 
meanings of the new parameters and their effects on accelerator performance were explored through a series 
of theoretical arguments and numerical simulations. The value of current rise rate through the coil is used 
as a proxy measure for the voltage across the coil, with higher values of dl/dt potentially yielding better 
inductive propellant ionization. The following insights were gained in this study. 

• As in the single capacitor case, there exists a value of the dynamic impedance parameter a for which 
thrust efficiency is maximized. This optimum corresponds to a matching of the driving circuits natural 
oscillation time scale to the residence time scale of the current sheet in the acceleration zone. 

• The efficiency and maximum dl^/dt* through the coil suffer for lower values of capacitance ratio C, 
with much of the charge transferring back and forth from C\ to C2 instead of going through the coil 
and performing work on the propellant. 

• The efficiency of the two-capacitor system approaches the theoretical efficiency of the single capacitor 
topology for C 1, with the commensurate value of dl^/dt* asymptoting to a value twice as large as 
that possible in the single capacitor case. 

• As C is increased, the gross behavior of the current waveforms for I* and match. The differences in 
these waveforms consists of a high-frequency modulation overlayed on the gross waveform. For a given 
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C the modulation is a function of n, which represents the position of C2 relative to C\ and the coil. 
The value of n greatly affects the level of modulation on I* but only slightly affects , implying very 
little change in the acceleration process as a function of n. 
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Introduction 



• Typically pulsed inductive thrusters limited to a simple RLC ringing circuit 

• A single data point (from 1971) indictated adding a second capacitor in the 
system downstream of the main bank offered potential performance 
improvement and a greater inductive voltage drop across the coil 

• Implies higher dl/dt through the coil 

• Greater dl/dt through the coil could result in more efficient and/or complete 
inductive breakdown/ionization of the propellant 

• In this study 

• Model the system with a second capacitor, nondimensionalizing to 
explore the global properties 

• Calculate the effect of the second capacitor on performance 

• Calculate dl/dt through the coil, using it’s value as a proxy measurement 
for the ability to inductively ionize the propellant 
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Circuit Model 



• Two current loops, coupled through the inductive acceleration coil 

L 0 



• Equivalent circuit 
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Circuit Model w/Additional Capacitor 



• Added new capacitor C 2 

• Stray inductance Z_ 0 , external resistance R e split between loops for / 1 and / 2 

• Constant n associated with how close C 2 is to C 1 
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Circuit Equations 



► Governing equations for current (from Kirchoff’s law) 

dh V1-V2- nReh 


dt 
dJ 2 
dt 

dh 

dt 


nL 0 


L C V 2 + ( L C I 3 + Mh) ^ - (1 - n)R e L c h - RyMld / {L c [(1 -n)L 0 + L c ] - M 2 } , 


!nt * 

Lc 


• Governing equations for voltage on the capacitors (from charge conservation) 

dV 1 __ h dV 2 _ h-h 

dt ~ ~ Cx dt C 2 ’ 


• Governing equation for the mutual inductance (empirical model) 

dM Lc f z \ dz 

dt 2zq eX ^ \ 2zq ) dt 
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Acceleration Model Equations 


Mass Accumulation Equation (snowplow model) 

m ( t ) — mo + / pav z dt 

Jt = o 


Momentum Equation 


Lcjl 

2 zq 


‘Snowplowing’ 

exp (- 2 / 2 0 ) = Pav 2 z + rn (t) 



dt 


EM Force 


Acceleration 


Non-dimensionalization 



• Initial non-dimensional terms 


it 

n 



T* 

i 2 


t* = 



y/LoC 


• Additional terms 


z 


* 


z 

ZO 



M* 


M 

Lc' 


Vl 

Vo’ 



^2 

Vo’ 


5lc 

v z = 




i h 


Z () 


mo = mo/mbit p* — poW m bit 
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Non-dimensional equations 



dl r 
dt* 
dll 
dt* 
d ll 
dt* 
dV{ 
dt* 
dVj 
dt* 
dA£ 
dt* 
dz* 
dt* 

dv* z 

dt* 

dm* 

dt* 


- (VC - vo - iu>i 

n 


dM 


L*V 2 * + + /£) — I£M*L*ip 2 - (1 - n) I^L*tpi 


. r , dIS „ dM < 

A/*— - + /o— — 
dt* 2 dt* 


dt 


- /££> 2 


/ ((1 - n) L* + 1) - (M*) ; 


A* 


= c(/r-/ 2 *) 


2 exp “t 




= v. 


a (A ) 2 exp (—z*) 


P*/( 2 *)K ) 2 



* 


p 7 (*> 


* 

2 
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Initial Conditions / Scaling Parameters 



Initial conditions 

n (o) - o, 

M* (0) = 1, 

IS (0) = 0, 

2 * (0) = 0, 

IS (0) = o, 

< (0) = 0, 

V? (0) = 1, 

m* (0) = -^2_. 

V 2 * (0) = 0, 

mbit 



Capacitance Ratio 


Dynamic Impedance 


Results 



Increasing & 
decreases swing in 
Ai*, increases ringing 
frequency 




D) 


1 . 5 - 
1 . 0 - 
0 . 5 - 
0.0 
- 0.5 - 
- 1 . 0 - 


:WWWW 


G) 


i — i — i — i — i — r 


1 . 5 - 

1 . 0 - 

0 . 5 - 


- 0 . 5 - 

- 1 . 0 - 


3 - 


i — r 

0 1 2 


3 4 5 

e 






n = 0 . 2 , C = 10 


i — i — i — i — r 


Increasing g increases 
dridt* across the coil 


Increasing n decreases amplitude and ringing frequency of I* and dl 2 */dt*, but does not affect the average 
waveforms 
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Performance as a function of & and a 




C)1000 



0.1 1 10 100 


• For fixed efficiency has a maximum 
as a function of a 

• Efficiency at a constant a increases with 
increasing asymptoting for ^ between 
10 and 100 

• Current rise rate increases with 
increasing asymptoting to a value of 2 
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Comparing performance to single-capacitor case 




• As ^approaches 100, efficiency 
approaches that of the single-capacitor 
case 

• Current rise rate increases with 
increasing asymptoting to a value of 2, 
which is 2x the fixed value for the single 
capacitor case 

• Performance in unaffected, while 
potential to ionize is enhanced 
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Conclusions 



• As in the single-capacitor case, there is an optimum dynamic impedance 
ratio a that maximizes performance 

• Efficiency, maximum dl 2 */dt* suffer at lower values of 0, with charge 
oscillating between C 1 and C 2 instead of passing through the coil 

• Efficiency of the two-capacitor system approaches the single-capacitor 
case for &» 1 , with the commensurate dl 2 */dt* growing to twice the value 
of dl*/dt* through the coil in the single-capacitor case 

• As C increases, the gross behavior of I* and \ 2 match. The differences 
consist of differing high-frequency modulations overlaid upon the gross 
waveforms. 
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